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We study a lattice realization of a Luttinger liquid interacting with a bath of quantum spins in 
terms of an antiferromagnetic S=l/2 Heisenberg chain, where each spin is also coupled to a a = 1/2 
Kondo spin degree of freedom. This model describes the low-energy spin dynamics in quasi one 
dimensional materials, where the electronic spins of the magnetic ions interact with those of impu- 
rities, nuclei and possibly other spin species present in their environment. For large ferromagnetic 
and antiferromagnetic Kondo interaction J' there are two phases corresponding to an effective spin- 
1 Heisenberg chain and a dimerised spin- 1/2 ladder, respectively. For weak Kondo interaction we 
establish that the Kondo interaction drives the system to a strong coupling regime. This suggests 
that J' = is the only critical point in the system. 

PACS numbers: 71.27.-|-a, 75.10.Pq, 75.10.Jm, 75.40.Gb, 75.50.Ee 



I. INTRODUCTION 

Entanglement and cooperative spin behaviour induced 
by coupling two spin systems with different intrinsic dy- 
namics is a recurrent theme in condensed matter physics. 
In metals, it emerges in the context of the coupling 
between the spins of itinerant electrons and localized 
atomic/impurity spins, known as Anderson impurity or 
Kondo problem^. The famous "central spin" problem of 
an electronic system interacting with a bath of nuclear 
or impurity spins has recently re-emerged in the area of 
spin-polarized microelectronics and the physics of quan- 
tum, spin-entangled electronic states^. Perhaps the most 
interesting example of this kind is a quantum-critical 
spin system, such as found in strongly correlated mag- 
netic insulators, coupled to a bath of quantum or classi- 
cal spins^iiii^ii^i^. In practice, the interaction between 
electronic spins in the system of interest and "external" 
spin variables arises in a variety of different contexts and 
occurs on vastly different energy scales. The first exam- 
ple is the exchange interaction between valence electrons 
involved in cohesion or chemical bonding and a localized 
"d" -orbital, such as in the s-d model of magnetic metals. 
If the d-orbital is singly occupied and the hybridization is 
weak, the dominant interaction is Kondo exchange^. Sec- 
ond is the coupling of electronic spins to impurity spins 
present in real material. While its effects scale with im- 
purity concentration, in many cases and in particular at 
low temperatures it is the determining mechanism for 
physical phenomena such as damping and quantum de- 
coherence. The response of a spin system that is either 
at or close to quantum criticality to impurity spins is 
particularly singular. It exhibits fascinating impurity- 
driven physics,-'—, which has emerged in studies of lightly 
doped cuprates and related two-dimensional Mott insu- 
lators. Another example occurs in complex alloys, where 



in addition to magnetic Sd ions there often exist mag- 
netic rare-earth cations (R'^"''), which lead to a lattice 
of macroscopically many "impurity" spins even for ideal 
stoichiometric materials^. In the Haldane (S=l) chain 
atiferromagnet R2BaNi05, the cooperative coupling to 
paramagnetic rare-earth spins dramatically modifies the 
spin dynamics of gapped Ni^"*" chains and induces mag- 
netic order at a finite temperature^. Finally, at very 
low energies/temperatures the hyperfine coupling of elec- 
tronic and nuclear spins becomes importantii^i^. Indeed, 
many abundant isotopes of magnetic ions have non-zero 
nuclear spin^''. Furthermore their electronic spins can 
also interact with nuclei of surrounding ligand ions^^. 

Although in many cases the coupling of the spin sys- 
tem to external spin degrees of freedom can be neglected 
in the same way as the coupling to a generic thermostat 
is swept under the carpet in equilibrium statistical me- 
chanics, it is of the same fundamental importance. The 
existence of such a coupling is required in order for the 
quantum (spin) system to reach its equilibrium state, or 
the ground state at T=0. Relaxation to equilibrium re- 
quires the change of both energy and angular momen- 
tum, which are integrals of motion for an isolated spin 
system. The simplest example is a Heisenberg magnet 
(or paramagnet) in a magnetic field. Here the Hamilto- 
nian conserves the total spin component along the field 
direction. It is only through (the implicit) coupling to 
some external system of angular momenta, or spin bath, 
that the magnet can adjust its total spin as the magnetic 
field changes (e.g., in a quantum phase transition from a 
spin-gap to a magnetized phase). 

While the effect of external spin degrees of freedom on 
the spin dynamics can often be averaged out, e.g., in the 
framework of spin-boson or spin-bath model o^^i^^ , there 
are important cases where the coupled dynamics of the 
two systems is of paramount interest. One such case is 
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known as "pulling" - which refers to the hybrid dynamics 
in a system of electronic spins coupled to a thermalized 
ensemble of nuclear spins. Interest in this phenomenon 
was recently renewed by studies of field-induced quan- 
tum phase transitions in quantum magnets, where it was 
discovered that pulling prevents the expected full soft- 
ening of the electronic spin excitation spectrum^^. In- 
stead, the latter acquires a gap, which increases with 
decreasing temperature, while a soft-mode behaviour is 
induced in the nuclear spins, which "take over" the quan- 
tum criticality^!^. 

Here we focus on the opposite, much less studied and 
understood limit of a quantum-critical spin system cou- 
pled to a bath of quantum spins, whose dynamics is gov- 
erned solely by their Berry phases. We consider a Heisen- 
berg S=l/2 chain with antiferromagnetic exchange cou- 
pling J > 0, where each spin is coupled to an additional 
local Kondo spin cr = 1/2 by an exchange J', 

N N 

H = JY,S{^)■S{J) + fY,S{i)■cr{^). (1) 

(ij) i=l 

This model is sometimes referred to as incomplete ladder, 
or SU(2) Kondo-necklace model, Fig. [TJ The model ([T]) is 
a one-dimensional (ID) analog of the incomplete bilayer, 
which recently received much attention in the context of 
the 2D cuprates^*^. The fundamental difference between 
the ID and 2D Kondo necklaces is that in D=2 quantum 
criticality is achieved by tuning the coupling J', while in 
D=l the S=l/2 chain is a critical Luttinger liquid for 
J' = 0. In what follows we allow the Kondo-coupling 
J' in ([H to the local spins to be either ferromagnetic or 
antiferromagnetic. 

S(l) _ S(2) S(3) 

---0- -0- -0-- 



!j' Ij' !j' 




FIG. 1: One-dimensional Kondo necklace model: a "master" 
S = 1/2 Heisenberg chain with antiferromagnetic coupling J, 
where each spin interacts with an extra, Kondo spin degree 
of freedom a = 1/2, via coupling J'. 

The model ([!]) is closely related to the Kondo necklace 
model, 

N N 

Hknm = jJ2s^i)S^j) + Sy{i)Sy{j) + J' ^ S{z)-a{i), 

_ . . (2) 

which was introduced in Ref. [l^ as a simplified ver- 
sion of the Kondo lattice model in the ID case. The 



Kondo necklace model ([2]) has been studied by a vari- 
ety of methods such as Monte Carlo simulationsii, real- 
space renormalization group (RG) techniques'^, exact 
numerical diagonalization^^, density-matrix renormaliza- 
tion group (DMRG) computations^*''^^, bond operator 
mean field theory22i2^ and bosonization^iiS^. While all 
methods agree in the regime | J'| ^ J, it is still controver- 
sial whether a small Kondo coupling \J'\ <C J leads to the 
formation of a spin gap. While real-space RG and exact 
diagonalization of small systems suggest the existence of 
a critical value J'^ of the Kondo coupling strength, below 
which the spin gap vanishes, Monte-Carlo and DMRG 
computations indicate that in fact = 0. The doped 
case has recently been investigated in^^. In the two and 
three-dimensional cases it is well established that there 
is a quantum phase transition between a Neel ordered 
ground state at small J' > and a Kondo-singlet phase 

at J' > 7 ^22.23.27.28.29 

In the present work we study the model ([1]) using the 
flow equation method. In order to "close" the system of 
flow equations we employ a decoupling scheme for spin 
operator products, which is asymptotically valid in the 
limit J' ^ 0, and use essentially exact expressions for 
the two-spin correlation function in the "master" spin- 
1/2 chain, derived by field theory methoda^i^. The flow 
equations obtained upon decoupling are then solved nu- 
merically in the weak coupling limit. 



II. STRONG COUPLING LIMITS 

In order to understand general properties and phase 
diagram of the model, it is instructive to analyze the 
strong coupling limits J' ~* ±oo first. 

A. Ferromagnetic Kondo coupling J' — oo 

For strong ferromagnetic Kondo coupling an effective 
spin 1 is formed on each rung. The Hamiltonian describ- 
ing the interaction between these spins takes the form of 
an antiferromagnetic spin-1 Heisenberg chain 

HoS = ^Y.^rT,+i. (3) 

j 

Here T" are spin-1 operators. The antiferromagnetic 
spin-1 chain is known to display a spontaneously bro- 
ken Z2 'Si Z2 symmetry characterized by a non-zero string 
order Oftj-ing 7^ 0' where^ 

/ j+n-l \ 

Os^ing - ( exp(z. ^^")^"+« ) • (4) 

\ k=j+l I 

It is well known that excitations in the spin-1 Heisen- 
berg chain are desribed in terms of a triplet of gapped 
magnons. 
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B. Antiferromagnetic Kondo coupling J' — > cxo 



In this limit the ground state is that of decoupled sin- 
glet dimers. This is in fact the same ground state as for 
the regular two-leg ladder with Jrung 3> Jiog in the limit 



J, 



rung 



oo . It also can be characterized by a string order 



parameter, e.g. 
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Ij+n 



The excitation spectrum in this limit again has a gap. 



(5) 



fe 



(9) 



Here the Fourier transformed spin operators are defined 
as 

Sfe-5^S(n)e-^^", S(n) = l^Sfce^'=". (10) 



III. WEAK KONDO COUPLING 

The question we want to address is what happens for 
weak Kondo couplings \J'\ <C J. The Kondo necklace 
model can be viewed as a particular limit of an asym- 
metric two-leg ladder model, in which the coupling along 
the first leg J is much larger than the rung coupling J', 
which in turn is large compared to the exchange J2 along 
the second leg 



J > J' > J2 



(6) 



This case is difficult to analyze for the following reason. 
Bosonizing the spin chains making up the two legs of the 
ladder results in a two-flavour Luttingcr liquid. However, 
the cutoff of this theory is equal to J2. The rung cou- 
pling J' can then not be treated as a perturbation of the 
two-flavour Luttinger liquid as it is much larger than the 
cutoff of the latter. In the Kondo necklace model the role 
of J2 is played by the RKKY interaction induced by J', 
which is of order J'^/J < | J'|. 

In order to analyze the small J' regime we have em- 
ployed Wegner's flow equation method2ii2^i^i22i^i2S. 



A. Flow Equation Method 

In the flow equation method^i^ a one-parameter fam- 
ily of unitarily equivalent Hamiltonians H(l) is con- 
structed via the differential equation 



dH{l) 
dl 



The anti-hermitian generator rj{l) is taken as 
r^{l) = [H^{l),H{l% 



(7) 



(8) 



where H^il) is a particularly chosen "diagonal" part of 
the Hamiltonian. For the Kondo necklace model we chose 
Ho{l) as 



In the initial Hamiltonian ^ the second term is ab- 
sent, but it will be generated under the flow. The full 
Hamiltonian will be of the form 



Hil) = Hoil) + Hi{l) + H2{l) , (11) 



where 



^iW = ^E^fcWSfc-'^-^- ■ (12) 



The contribution H2{1) will loosely speaking contain all 
multi-spin interaction terms compatible with the global 
SU(2) spin rotational symmetry. A key element in im- 
plementing the flow equation approach is that the initial 
coupling J' is small and all terms in i?2(0 "^iH be of order 
J'^ or higher. As long as we constrain our attention to 
the small J' limit, we may therefore neglect i?2(0 when 
calculating the generator r]{l) of the unitary transforma- 
tion 



r?(0-[i?o(0,^fi(0]=m(0 + ^2(0- 



The explicit forms of ryi^2(0 ^-^^ 



(13) 



Viil) 



mil) 



i 



kk' 



[S 



^4,(0K(0-«fe+fe'(0] 



kk' 



(14) 



Working out the required commutators of f?i,2(0 with 
-ffo,i(0 (see Appendix E]) we arrive at the following ex- 



pression for the Hamiltonian H{1) 
H{1) = Jk{l)Sk ■ S-fc + ^ ak{l)eTk ■ c-k 

k k 
-^^J'k{l)^k ■ (T-k 



J3 E<P..W Sfc+p+,-^-p S. 



iV3 



k,p,q 



J2 ^^kl,q{^) Sfc+p • S_fc • cr_p_ 



k,p,q 



iV3 



k,p,q 



k,p,q 



p-<? 



k,p,q 



J2 ^''^kl,qi^) Wk+p+q X (T_fe] • S_p_5. 



k,p,q 



(15) 



The initial values of the various couplings are 

Jfe(O) = Jcosfc, 4(0) = J', Mj^l^{0) = 0. 

In order to obtain a set of flow equations we need to de- 
couple the three- and four-spin terms. We do this by ex- 
panding in fluctuations around the J' = ground state, 
i.e., 

Sk ■ Sk' (Tq ■ (T-k-k'-q — k ' ^ k') (J' q ' -k-k' -q 
+ {cTq ■ (T^k-k'-q)Sk ' Sfc' 
+ : Sfc • Sfc' (Tq ■ CT-k-k'-q ■ , (16) 



Sk' ■ (Tq Sk ■ CT-k-k'-q — g (Sfc' • Sfc)(Tg • Cr_fe_fc^_q 
+ 7:{'^q • 0--k-k'-q)^k' ■ Sfc 



Sfc ■ Sk' Sq ■ (T^k-k'-q 



(Sfc • Sfe')Sg • (T-k-k'-q 



--(Sfc • Sq)SA;' • Cr^k-k'-q 
"^(S/c' • Sq)Sk ■ (T^k-k'-q 

- : Sfc • Sfe' Sq • cT_i._i.i_q : . (18) 



These expansions can be motivated by bosonizing the 
Hamiltonian (see Appendix |B| . Substituting (flS)) 
into we obtain expressions for Hq{1) and Hiil). The 
static spin-spin correlation functions entering (fT5|) and 
hence the expressions for ifo,i(0 should be calculated 
self-consistently with respect to the flowing Hamiltonian 
Ho{l). However, at weak coupling | J'| ^ J one may cal- 
culate the correlators with respect to the initial Hamilto- 
nian Hq(0), as the corrections are of higher order in J'. 
Taking into account that 



(Sfc • Sfc') — 5fc^_fc' (Sfc • S_fc) 



(19) 



and retaining only terms quadratic in spin operators, the 
Hamiltonian ITSl) is reduced to 



^ -^Y'^kH) ^k ■ ^-k + ^^Oik{l) (Tk ■ (T-k 
k k 

+ ^J2-^kil)Sk-rT_k. (20) 



Sfc' • (Tq Sfc • (T^k-k'-q '■ , (17) Thc flow equatious for the couplings take the form 



dJk _ 2 \^(p.2fr,j J J ^ {f^k- ■ cT-k') , ^J'kf 



/O ■ ^-k') '0-1 N 

(2afc' - Qifc+fc' - afc-fc') , (21) 



dak 2 



dl 



3N 



Y^iJ'yfi-^ak - ak+k' - ak-k') ^^""^-"'^ + ^^T.(^Jk' ~ Jk+k' ~ Jk-k' 



N 



(22) 



5 



= '^'^{(Jk' - Jk){Jk+k' + Jk-k' -'2Jk , 
k' 

+ ^ 2^'Jk'Jk+k'\[Jk+k' - Jk) 



- + [ak' - a/c ) (ak+k' + ctk-k' - i-aw ) 



N 



+ {ctk+k' 
I 



afe) 



\C k+k' ■ CT-k-k 



N 



(23) 



As we have remarked earlier, in order to solve the flow 
equations we need to know the static spin-spin correla- 
tor {Sk ■ S_fe) for the isotropic spin-1/2 Heisenberg chain. 
This can be calculated accurately from the results of Refs. 
"30ll3l| . There the large distance asymptotics of spin- 



spin correlations functions was determined by combin- 
ing exact results for correlation amplitudes^ with renor- 
malization group improved perturbation theory in the 
marginally irrelevant interaction of spin currents present 
in the continuum description of the spin-1/2 Heisenberg 
chain that 



{Sim + j)-S{j)) 



where 



(-1) 



h{g) 



(24) 



5 
128 



c 

16 



f2{g) = i + f + (c + 



21 
1024 
9_ 
2 



256 





c3 


13C(3) 


64 


16 + 


32 


^)^ + 


c{c + 




4/2 


2 





Here g is the running coupling constant 



(25) 



(26) 



and c is a free parameter that related to the choice of 
renormalization scheme. It was demonstrated in Ref. [3l| 
that for I > 2 p4|) is in very good agreement with numer- 
ical results. Supplementing these results with the known 
values 



(S(j)-S(j + 1)> 
(S(j)-S(j)) 



n=l 



-0.443147.. 



n 



(27) 



and then Fourier transforming, we arrive at the result 
shown in FigO We see that in momentum space the 
spin-spin correlator is dominated by the logarithmic di- 
vergence at = TT. 



B. Solution of the Flow Equations 

We are now in a position to solve the flow equations 
(PT|) - (|23p numerically. We find that for small initial val- 
ues on the Kondo interaction \ J'\ < O.IJ the behaviour of 
the solutions to the flow equations exhibits two regimes. 




FIG. 2: Static spin-spin correlation function on the isotropic 
spin-1/2 Heisenberg chain. 



In the following we measure I in units of J~^. For large 
flow parameters I > 100 the couplings appear to ap- 
proach a fixed point. However, as I increases further, 
the behaviour eventually changes and the couplings di- 
verge. We interpret the eventual runaway flow as being 
indicative of the emergence of a strong coupling phase, 
characterized by the formation of a spectral gap. 



1. Kondo Interaction 

We find that for I > 100 and small initial values of the 
Kondo interaction | J'| < O.IJ, the couplings become 
very small everywhere except at k ~ 0, 2tt. This holds in 
the ferromagnetic as well as the antiferromagnetic case. 
This behaviour is shown in Figl3l where we plot as 
a function of k for different values of I and J' — 0.05J. 
The vicinity of fc = 27r is shown in greater detail in FigH) 
In Fig Owe show the Kondo coupling in the vicinity of 
k = n for several values of the flow parameter I. We ob- 
serve that for large values of I J'i.{l) appears to approach 
a small but finite limit. In position space the Kondo cou- 
pling J' becomes more and more long-ranged as the flow 
parameter I increases. This is shown in FiglH] For very 
large values of the flow parameter I we always enter a 
regime in which the diverge. This is shown for the 
case J' = 0.06 J in FigH 
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FIG. 3: Kondo coupling J'f. as a function of the momentum 
k for several values of the parameter / characterizing the flow 
; = 1, ; = 4, / = 8, Z = 32 and ; = 128. The initial value is 
Jj. = 0.05 J. The Kondo interaction becomes small under the 
flow except in the vicinity of A; = 0, 27r. 



FIG. 5: Kondo coupling as a function of the momentum k 
in the vicinity of A; = tt for several values of the parameter I 
characterizing the flow. The initial value is — 0.05J. The 
couplings J'f. tend to a small but finite limit for large values 
of L 



0.05 



0.02 





FIG. 4: Kondo coupling J'f. in the vicinity of A; = 27r for several 
values of the parameter I characterizing the flow. The initial 
value is J'k = 0.05J. 



FIG. 6: Kondo coupling J'(n) in coordinate space for valious 
values of the flow parameter The initial value of the Kondo 
coupling is J' = 0.05J . There is a long range smooth com- 
ponent as well as a shorter range staggered one. 



2. RKKY Interaction 



3. Heisenberg Interaction 



The flow of the RKKY couplings ak{l) is shown in 
Figl5]for several values of the flow parameter I. We see 
that the flow appears to approach a small flnite limit 
quite quickly. The RKKY couplings are strongest in the 
vicinity of fc = vr, where they become of order J'^ / J. 

For sufRciently large values of the flow parameter / the 
RKKY interaction starts to decrease quite quickly and 
eventually diverges. 



The Heisenberg interaction changes only weakly under 
the flow. We therefore plot the difference Jk{l) — Jk{0) 
rather than Jk{l) itself in FigO We find that the k de- 
pendence is unchanged under the flow, so that 



Jk{l)^J{l)cos{k) 
Here J{1) ^ J + 0{ j'^ / j) . 



(28) 
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FIG. 7: Kondo coupling as a function of the momentum 
k in the vicinity of A; = vr for large values of the parameter I 
characterizing the flow. The initial value is J' = 0.06J. The 
absolute values of the couplings are seen to increase with I. 



FIG. 9: Difference of Heisenberg couplings Jk{l) — Jk{0) as a 
function of the momentum k for several values of the param- 
eter I characterizing the flow I — 1, I = 4, I — 8, I = 32. The 
initial value was taken to be — 0.05J. 
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FIG. 10: Zero temperature phase diagram of the SU(2) in- 
variant one dimensional Kondo necklace model. 



does not allow for a determination of how A scales with 
J'. 
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FIG. 8: RKKY coupling at as a function of the momentum k 
for several values of the parameter I characterizing the flow. 
The initial value is = 0.05J. 



IV. DISCUSSION 

The above analysis of the flow equations shows that 
a small Kondo interaction drives the system to a strong 
coupling regime. On the basis of the analysis presented 
in this work we cannot establish the nature of the strong 
coupling phase. However, using exact diagonalization of 
small systems one can see that there is a spin gap for J' > 
0.3J, which suggests that the large- J' phase extends at 
least down to such interaction strengths. Combining this 
observation with the result of the flow equation analysis 
suggests a form of the phase diagram as shown in FigfTUl 
This would imply there presence of a spin gap A for any 
J' 0, but the flow equation approach presented here 
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APPENDIX A: COMMUTATORS 

In this appendix we list the commutators needed for 
determining the flow equations. 

[m,Ho] ^ j;pYl ^S,<2 ^ ^-'^ ' ^ '^-9 

k,p,q 

+ ]^ E ri'Ls, • s_fc Sfe+,_, . (Ai) 

k,p,q 
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APPENDIX B: BOSONIZATION 



\n fT ] — ^ \^ r^'^) ^, V , c! V c; The commutator [772,^11 is obtained from [771, iJi] by re- 

" l^^f'^P^'^'''^''"'-'-^'-'''^-'^ placing r(5) by r(6) and interchanging 5" I; a". 

+ E ■ '^-^ '^k+p-q ■ S-p. (A2) 

k,p,q 

— ^k+p+q • (T_g • Cr^p 

-i [Sfc+p+g X S_fc] • cr_p_q + h.c.|.(A3) 

In order to bosonize the Hamiltonian (jlSp it is con- 

^ venient to tranform the spin operators to coordinate 

^k,p,q ~ Jp[Jk — Jk+p][o:q — ap-q], space. At weak couphng the generated interactions are 

p(2) _ —4 J' \Ji^ — Jj. 1 [J — Jj. ] short-ranged. In a continuum description the most rele- 

'^'P'? P '^^ ^ vant (in the renormahzation group sense) contributions 

v'f'^ — Jp [ak — ctk+p] [Jq — Jp~q], can then be obtained by considering nearest-neighbour 

f^^j ^ terms. Interactions between next nearest neighbours will 

^k,p,q ~ ~4Jp [ak — ctk+p] [Jq — Jk+p-q], merely lead to a renormahzation of the couplings as well 

p(5) _ —J' J' [Jk — Jk+ ] generate (less relevant) derivative terms. Using PH)) 

'^'^''^ ^ ' p ' keeping only the terms corresponding to nearest- 

^'kp q ~ ~J'qJp V^k — oik+p] ■ (A4) ucighbour interactions we find 



H{1) - J{l)^S{i)-S{i^-\) + a{l)^cT{i)-a{i + l) 

i i 

+ J'(/) ^ S(z) • <T(») + J' {I) S(z) • {cT{i + 1) + a{i 1)) 

i i 

+ M(i)(0 ^ S{i) ■ CT{i + 1) S(^ + 1) • c^(^) + M^^) (?) ^ S(i) • S(i + 1) S(^ - 1) • (cr(i) - cr(i + 1)) 

z i 

+ M^^^il) J2 S(«) • S(^ + 1) cr(z) • cr(i 1) -|- M^-*) (0 ^ cr(z) • cr(i + 1) cr(i - 1) • (S(^) - S(i + 1)) 

i i 

+ zM(5)(/) ^[S(i) X S{i + 1)] • (*T(i) - (t{i + 1)) + zM(6)(/) ^[cr(i) X ,T(^ + 1)] • (S(^) - S(^ + 1)). (Bl) 



I 

We now bosonize the part = J it) S(i) • S(i + Here z = ur — ia;, z = vt + ix, ^ = Lp + (p, O ~ ip — (f 

1) of the Hamiltonian (jBip by standard methods, see and we use a normalization such that 

e g ^33.41.42.43.44 

5"(j) J"(z) + J^-z) + (-l)%"(x) , (B2) The bosonized form of i/o is 

^ ^e--W , J+(z-) = ^e'^(^) , ^" = lB f + 

V ^ 27r ^ ^ 27r -^O'^ 

J^{z) = —i—dzf , •^^(^) ~ "^^2"? ' Using operator product expansions we then can extract 

Q Q ^ the dominant parts of the various other terms in (|Bip . 

n(a;) — c-yao(^cos(— ), — sin(— ), — sin(— ) j .(B3) Denoting x = jao we have 
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S(j).S(j + 1)<t(j).^(j + 1) 
S(j).^(j + 1) S(j + l).^(j) 



ci + (-!)■' C2 cos(^) cr(j) • cr(j + 1) + . . . , 

ci + (-1FC2 cos(|)]cr(j) • cr(j + 1) + C3[J(.t) - J(a;)] • cr(j) 



(S(j) X S(j + 1)) • <t(j) ^ 2C3 [J(x) - J(x)] • aij) 



(B6) 

(B7) 
(B8) 



S(j) -SU + l)S{j - 1) • - S{j - 1) -Simj - 2) • <t(j) cx (-l)^n(a;) • ^(j). 



(B9) 



As expected, the interactions generated under the flow 
are simply all terms that are compatible with the global 
spin rotational SU{2) symmetry. The term (J(x) — J(a;))- 
cr(j) breaks the left-right (chiral) symmetry of the free 
boson Hamiltonian in the J' — > limit. However, it 
is well known^ that interchange of left and right mov- 
ing bosons is not a symmetry of the spin-1/2 Heisenberg 
model at low energies due to the presence of a marginally 
irrelevant interaction of spin currents (which we have 
omitted in (|B5|1 for the sake of brevity). Inspection of 
the scaling dimensions of the bosonic parts of the vari- 
ous terms generated along the flow suggests that at low 
energies the most relevant interaction of staggered mag- 
netizations is precisely the one picked out by our decou- 
pling scheme (IT51) . We note that the term (|B8p is not 



included in our decoupling scheme. While it is likely to 
be as relevant as J" ■ <t, the bare coupling of the latter is 
much larger, which justifies neglecting the former in the 
weak coupling regime J' <C J. 

Another fluctuation induced contribution that is not 
included in our decoupling scheme is the interation of 
dimerizations 

{-iycos{^) a{j)-rT{j + l). (BIO) 

It is less relevant than the terms we keep. If it were the 
dominant interaction in the problem, its effect would be 
to open up dimerization gaps among the Kondo spins 
and spin-chain spins respectively. 
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